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IV 


1  SUMMARY 


The  University  of  Missouri  at  Columbia,  MO  and  the  Air  Force  Research  Laboratory, 
Kirtland  AFB  NM,  collaborated  on  this  research  study  pertaining  to  the  fabrication,  analysis  of 
scaled  down  semiconductor  devices,  and  theoretical  quantum  mechanical  modeling.  The  work  was 
carried  out  under  the  AFRL  Grant  No.  FA9453- 10- 1-0001.  The  semiconductor  processing, 
fabrication  and  the  resulting  carrier  transport  characteristics  of  Mobile  Service  Manager  (MSM) 
devices  fabricated  as  wall  like  structures  in  silicon  on  insulator  technology  were  reported.  MSM 
device  dark  current,  direct  current  (DC)  photocurrents,  and  the  time  response  of  carrier  transport 
were  investigated.  The  resulting  conducting  channels  were  actually  smaller  than  their  physical 
dimensions,  a  result  of  depletion  of  carrier  near  the  interfaces.  As  the  physical  channel  widths  were 
reduced  by  oxidation,  strain  was  produced  near  the  interface  and  strained  lattice  became  a 
significant  portion  of  the  conducting  channel.  The  increase  in  mobilities  for  both  holes  and 
electrons  stemming  from  the  strained  silicon  resulted  in  a  dramatic  increase  in  carrier  mobility  for 
both  electrons  and  holes  as  the  physical  channel  width  was  reduced  from  200  nm  to  20  nm.  The 
theoretical  model  incorporating  the  effects  of  strain  present  in  these  nanoscale  MSM  devices 
compared  favorably  with  experimental  results,  showing  that  hole  mobilities  increased  with 
decreasing  L.  Additionally,  if  these  electron  and  hole  mobilities  can  be  retained  with  the 
application  of  gate  electrodes,  then  this  technique  may  yield  a  much  simpler  path  towards  high 
performance  Complementary  Metal  Oxide  Semiconductor  (CMOS),  both  n-channel  and  p- 
channel,  than  current  techniques  for  either  planer  ultra -thin  body  Field  Effect  Transistor  (FET)s  or 
Fin  shaped  Field  Effect  Transistor  (FinFET)s. 
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2  INTRODUCTION 


In  thin  film  silicon-on-insulator  (SOI)  technology  a  buried  oxide  layer  separates  the  active 
Si  device  structure  from  the  parasitic  influence  of  the  substrate,  thus  making  the  device  less 
susceptible  to  radiation  effects,  specifically  single  event  upset.  [1]  The  fabrication  of  such  devices 
is  fully  compatible  with  existing  very-large-scale  integration  (VLSI)  fabrication  technology.  [2] 
Moreover,  it  has  become  possible  to  develop  and  study  electrically  isolated  ultra-thin-films  and 
nanostructures  for  the  quantum  confinement  effects  on  electrons  and  other  quasi-particles,  such  as 
holes  and  excitons.  [3] 

The  effects  of  quantum  confinement  on  carrier  transport  properties,  however,  have  been 
primarily  investigated  in  ternary  and  quaternary  material  heterostructures  and  superlattices,  in 
which  scattering  is  seen  to  enhance  some  modes  of  electron-lattice  interaction  while  suppressing 
others,  thereby  changing  the  relative  value  of  the  carrier’s  effective  masses  of  electrons  and  holes, 
as  compared  to  bulk  semiconductors.  [4]  Such  studies  in  Si  have  been  very  limited. 

In  the  nano-scale-regime,  electrons  start  to  behave  more  like  waves  than  particles,  yet  this 
transition  to  quantum  mechanical  regime  does  not  come  about  abruptly.  Rather,  there  is  a  transition 
region  in  which  bulk  properties  begin  to  slowly  weaken  while  the  quantum  mechanical  effects 
begin  to  strengthen.  [5]  Because  of  the  lack  of  available  scientific  literature  on  this  subject  in  Si,  it 
is  evident  that  little  attention  has  been  focused  in  this  transition  bulk/quantum  region  of  Si. 
Nevertheless,  an  understanding  of  the  physics,  in  particular  electron-hole  transport  properties, 
including  the  role  of  biaxial  strain,  in  the  transition  region  of  the  semiconductor  devices  is  of  great 
importance  due  to  the  increased  interest  as  the  VLSI  devices  are  scaled  down  to  the  nanometer 
regime.  [6] 

Strain  is  an  important  aspect  of  CMOS  and  has  been  incorporated  into  planar  CMOS  as 
standard  practice  since  the  late  1990s,  specifically  in  the  90  mn  technology  node,  in  order  to 
enhance  carrier  mobilities  and  ultimately  device  drive  currents.  As  detailed  by  S.  W.  Bedell  et.  al. 
[7],  biaxial  wafer  scale  tensile  strain  was  accomplished  first,  achieved  by  growing  silicon 
epitaxially  on  relaxed  GexSii-x  since  the  lattice  constant  of  silicon  was  less  than  the  germanium- 
silicon  alloy.  While  this  silicon  germanium  wafer  scale  approach  resulted  in  enhanced  nFET  drive 
current,  there  was  no  discernable  p-type  field-effect  transistor  (pFET)  drive  current  enhancement. 
By  incorporating  GexSii-x  layers  in  the  source  and  drain  regions  of  pFET  devices,  uniaxial 
compressive  strain  was  achieved  and  resulted  in  50%  enhancement  in  hole  mobilities  and 
subsequently  p-type  field-effect  transistor  (pFET)  drive  currents.  The  result  was  a  biaxial  tensile 
strained  n-type  field-effect  transistor  (nFET)  and  a  uniaxially  compressive  strained  pFET 
fabricated  simultaneously  in  the  CMOS  application. 
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Continued  scaling  of  CMOS  in  accordance  with  Moore’s  Law  and  the  International 
Technology  Roadmap  for  Semiconductors  (ITRS)  roadmap  required  increases  in  both  uniaxial  and 
biaxial  strain  to  achieve  necessary  drive  currents.  Uniaxial  strain  was  augmented  by  increasing  the 
Ge  content  in  the  source  and  drain  GexSii-x  layers,  closer  proximity  of  the  gennanium  bearing 
layers  to  the  channel,  and  with  new  gate  structures  designed  specifically  to  transfer  strain  to  the 
channel.  Uniaxial  tensile  strain  was  also  achieved  in  the  channels  of  nFET  devices  as  well  by 
incorporating  a  layered  source  and  drain  using  carbon  rather  than  germanium  to  reduce  the  lattice 
parameter  of  silicon. 

As  device  pitch  was  reduced,  localized  stressors  for  pFETs  became  less  effective  due  to  the 
reduced  volume  of  GexSii-x  material,  and  the  stressor  for  nFETs,  carbon  bearing  silicon, 
approached  the  limits  of  solubility  of  carbon  in  silicon.  Ultimately,  the  practical  limit  for  this  type 
of  strain  scaling  may  be  reached  as  a  result  of  severely  reduced  circuit  lifetimes  or  premature 
material  failure  [7]. 

A  scalable  approach  achieved  by  utilizing  wafer  level  strain  with  devices  fabricated  on 
mesas  providing  a  virtual  local  substrate,  has  proved  beneficial  for  device  pitch  below  100  nm. 
Strained  virtual  substrates  are  grown  on  insulator  material,  germanium-silicon  provides 
compressive  strain  for  pFETs  while  strained  silicon  on  insulator  provides  tensile  strain  for  nFETs. 
Initially  strain  is  biaxial  but  as  mesa  widths  are  reduced  below  0.3  pm  biaxial  strain  is  converted 
to  beneficial  uniaxial  strain.  While  this  provides  a  promising  path  towards  reduced  pitch,  edge 
relaxation  still  poses  a  problem  for  device  reliability  [7]. 

Regarding  three  dimensional  devices  such  as  new  multi-gate  or  tri-gate  architectures,  also 
known  as  FinFET  technology,  although  these  devices  have  achieved  performance  enhancement  by 
virtue  of  increased  sub-threshold  slope  and  a  reduction  in  drain  induced  barrier  leakage  with  the 
gate  geometries,  they  are  operated  at  relatively  high  loff  and  higher  voltage  drain  supply  (VDD) 
to  achieve  higher  switching  speed.  At  these  increased  operational  parameters,  higher  channel 
mobilities  becomes  desirable.  Process  induced  stressors  in  the  form  of  embedded  silicon 
gennanium  layers  in  the  source  and  drain  have  been  modeled  and  shown  to  have  a  modest 
performance  increase,  approximately  one-half  of  the  mobility  boost  for  similar  size  planar  FETs 
[7].  It  turns  out  that  process  induced  strain  in  FinFETs  is  most  effective  directly  under  the  gate  but 
diminishes  below  the  gate  and  along  the  sides  of  the  channel,  thus  reducing  average  mobility  and 
the  effectiveness  of  stressors.  Attempts  to  increase  strain  in  the  channel  and  along  the  sides  of 
FinFETs  with  the  choice  of  gate  material  and  oxides  are  limited  in  their  selection  by  the  work 
function,  inversion  capacitance,  fabrication  complexity  and  reliability.  Incorporation  of  wafer  level 
strain  using  SiGe-on-insulator  (SGOI)  and  Strained  Silicon-On-Insulator  (sSOI)  in  small  pitched 
circuits  may  be  possible  by  converting  the  tensile  strain  of  sSOI  to  compressive  strain  by  selective 
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growth  of  silicon-germanium.  The  industry  is  left  at  a  point  where  device  complexity  will  have  to 
be  increased  to  provide  methods  for  strain  incorporation. 

In  this  report  we  describe  the  comprehensive  study  where  the  nature  of  carrier  transport 
properties  are  evaluated  as  a  function  of  dimensional  scaling  wall  widths  from  200  nm  to  20  nm 
in  devices  where  the  width  is  reduced  by  thennal  oxidation  and  stress  naturally  accumulates  in  the 
device  channel.  We  believe  that  these  wall  structures  are  a  useful  starting  point  for  a  broader  study, 
as  these  can  be  configured  into  novel  high  density  3-D  devices,  where  thermal  effects,  such  as  heat 
buildup,  can  also  be  efficiently  managed.  [8] 
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3.  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

3.1  Rationale  for  substrate  material 

SOI  wafers  with  <  100  >  top  active  layer  crystal  orientation  were  used  to  fabricate  the  wall- 
like  structured  test  devices  for  this  study.  The  initial  SOI  structure  had  a  1500  rnn  top  active  layer 
on  a  3000  mn  buried  oxide.  The  SOI  configuration  allowed  complete  electrical  isolation  of  the  top 
active  Si  layer  from  the  underlying  substrate.  All  five  samples  in  this  study  had  identical  p-type 
active  layer  with  a  lightly  doped  concentration  of  1014  cm"3  boron  atoms.  Intrinsic  SOI  wafers 
would  have  been  an  ideal  choice  for  the  experiment,  however  due  to  the  commercial  unavailability 
of  completely  100%  intrinsic  material,  the  above  choice  of  dopant  type  and  concentration  was 
carefully  made  in  order  to  minimize  the  effects  of  impurity  scattering  as  a  result  of  redistribution 
of  impurities  during  the  thermal  oxidation  process.  Boron  was  chosen  as  dopant  because  it  tends 
to  segregate  away  from  the  Si  interface  and  into  the  thermally  grown  oxide,  thus  reducing  the 
impurity  concentration  near  the  Si  interface  with  SiCk.  [8]  Thus  the  segregation  coefficient,  which 
is  defined  as  the  ratio  of  the  dopant  concentrations  at  the  interface,  is  less  than  one.  This  process 
leads  to  the  formation  of  an  oxide  trapped  charge  (Qot),  which  contribute  to  the  formation  of  a 
depletion  region  near  the  Si/SiCk  interface.  [9,10] 

This  oxide  trapped  charge  (which  can  be  calculated  from  the  doping  concentration)  coupled 
with  the  estimate  of  fixed  charge  (Qf )  which  results  from  the  excess  Si  atoms  not  reacted  with  the 
oxygen  and  the  estimate  of  interface  trapped  charge  (Qit)  which  results  from  the  mismatch  between 
the  number  of  atomic  bonds  in  the  Si  crystal  surface  and  the  number  of  available  bonds  in  the  SiCk 
layer  combined  (Qot  +  Qf  +  Qit)  significantly  reduce  the  dimensional  cross-sectional  area  of  the 
channel  width.  10  Based  on  the  literature,  for  the  given  parameters,  it  is  reasonable  to  expect  that 
the  depletion  region  in  Si  would  extend  to  a  thickness  of  approximately  40  mn  from  the  Si/SiCk 
interfaces,  considering  the  thermally  grown  oxide  was  at  least  50  mn  thick.  [11]  Thus  a  physical 
Si  wall  width  of  approximately  100  mn  would  have  an  effective  width  of  20  mn. 

3.2  Crystalline  silicon  wall  nanostructure  fabrication 

In  order  to  fabricate  the  wall  structures,  photoresist  nano-scale  pattering  was  required.  The 
precursors  to  the  wall  structures  were  patterned  using  interferometric  lithography  (IL)  and 
reactive-ion-etching  (RIE).  IL  is  a  well-developed  technique  for  inexpensive  nano-patteming 
process.  [12]  IL,  in  its  simplest  fonn,  is  interference  between  two  coherent  waves  resulting  in  a  1- 
D  periodic  pattern  defined  by  7J2  sin  9  where  A.  is  the  optical  wavelength  and  29  is  the  angle 
between  the  interfering  beams.  A  typical  IL  configuration  consists  of  a  collimated  laser  beam 
incident  on  a  Fresnel  mirror  (FM)  arrangement  mounted  on  a  rotation  stage  for  period  variation. 
[13]  There  is  no  z-dependence  to  an  IL  exposure  pattern,  which  is  limited  only  by  the  laser 
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coherence  length  and  beam  overlaps.  [14]  The  1-D  nano-scale  pat-tems  were  first  formed  in  the 
photoresist  followed  by  pattern  transfer  onto  the  underlying  substrate  using  RIE  in  a  parallel  plate 
reactor  using  SF6  plasma  chemistry.  Figure  1  shows  a  scanning  electron  microscope  (SEM)  cross- 
sectional  image  of  an  array  of  nano-wall  structures  with  a  remaining  layer  of  patterned  photoresist 
after  RIE  has  been  performed. 


Figure  1.  Scanning  electron  microscope  (SEM)  cross-sectional  image  of  an 
array  of  wall  precursor  structures  with  a  remaining  layer  of  patterned 
photo-resist  after  reactive  ion  etch  process  step 

Note  at  this  stage  these  structures  are  merely  the  precursors  to  the  thin  wall  structures  that 
are  then  reconfigured  into  MSM  devices.  After  the  photoresist  was  re -moved  the  wall  structures 
were  thermally  oxidized,  in  the  absence  of  Fb,  in  order  to  further  minimize  any  effects  of  dopant 
(boron)  redistribution.  The  oxidation  process  accomplished  two  things.  First  it  consumes  the  Si 
thus  thins  the  wall  width.  Secondly  the  thermally  grown  oxide  preserves  a  low  defect,  clean 
Si/Si02  interface  and  at  the  same  time  passivates  the  surfaces  of  the  nanostructures  [15,16].  While 
the  Si/SiCb  interface  is  low  defect,  it  is  important  to  note  that  strain  is  present  at  the  interface  and 
reduces  with  distance  from  the  interface.  This  reduction  in  strain  as  a  function  of  depth  has  been 
seen  experimentally  in  Si/SiCb  interfaces  using  a  scanning  transmission  electron  microscope  using 
Z-contrast  imaging  to  produce  strain  contrast  imaging  [17].  Using  this  technique,  the  1/e  decay 
length  was  measured  at  0.2  nm.  The  modeling  of  the  thermal  oxidation  parameters  needed  for  the 
desired  thicknesses  was  complicated  due  to  the  fact  that  in  a  three  dimensional  wall  structure  there 
are  several  crystal  lattice  orientations  that  have  different  thennal  oxidation  rates.  As  a  first  order 
approximation,  we  used  average  values  of  oxidation  rates  between  the  various  lattice  orientations 
i.e.  oxygen  flow  rate,  pressure,  temperature  and  time.  These  parameters  were  then  fine-tuned 
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empirically  during  the  actual  thermal  oxidation  runs.  Figure  2(a)-2(c)  show  SEM  images  of  the 
cross-sectional  views  of  the  wall  structures  after  the  respective  thermal  oxidations. 

As  can  be  seen  from  the  SEM  images,  due  to  the  high  aspect  ratio  of  these  structures  the 
oxidation  rate  was  not  fully  unifonn  throughout  the  height  of  the  walls.  The  rate  was  faster  at  the 
top  part  of  the  walls  and  slower  at  the  bottom  part  due  to  higher  availability  of  oxygen  atoms  in 
the  upper  regions.  The  wall-like  structures  then  formed  the  active  region  of  the  MSM  devices  as 
described  in  the  next  section. 


Figure  2.  SEM  cross-sectional  images  of  any  array  of  wall  structures  after  thermal  oxidation;  (a)  200nm 
wall  structures  (b)  95nm  wall  structures;  (c)  40nm  wall  structures 


3.3  MSM  device  fabrication 

The  wall  structured  samples  were  then  configured  into  two  tenninal  metal- Si/nano  wall- 
metal  (MSM)  devices  for  optical  and  electrical  characterization.  The  MSM  device  configuration 
was  specifically  designed  so  the  current  would  flow  alongside  within  the  wall  boundaries  between 
the  electrodes.  This  allowed  the  physical  cross-section  of  the  wall  structures  to  dictate  the  current 
flow  properties.  The  mesa  structures  were  fabricated  to  cutoff  any  stray  current  paths  that  could 
bypass  the  intended  active  region  (wall)  carrier  path.  Figure  3(a)  shows  a  SEM  picture  of  a  typical 

n  1 i  ii  i 1  iTi  uffTi 

Figure  3.  SEM  image  of  (a)  pre-oxidized  Si  mesa  configuration  with  precursors  to  wall  structures  in 
the  active  region  in-between  planar  un-textured  regions  where  the  metal  contacts  will  be  deposited;  (b)  planar 
un-textured  Si  where  thermally  grown  oxide  was  removed  from  metal  contact  deposition  connection  walls;  (c) 

fully  fabricated  wall  device  with  interdigitated  electrodes. 
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pre-device  mesa  structure.  After  the  walls  were  oxidized  to  achieve  the  desired  wall  width,  the 
thermally  grown  oxide  was  selectively  removed  from  the  planar  un-textured  Si  pad  locations 
[Figure  3(b)]  using  an  appropriate  photo-mask  and  a  chemical  1:6  buffered  oxide  etch  (BOE) 
process. 

Following  the  resist  removal  the  samples  were  cleaned  using  a  sulfuric-acid:hydrogen- 
peroxide  solution  and  a  DI  water  rinse  followed  by  a  nitrogen  gas  dry  step.  The  samples  were  then 
re-patterned  using  photoresist  and  a  second  mask  was  used  for  the  process  to  form  the  electrode 
contact  regions.  Three  separate  evaporations  (30  mn  of  Ni)  were  performed.  The  first  one  was 
performed  at  a  normal  incidence  to  the  sample  surface  and  the  other  two  at  a  30o  degree  tilt  angles 
in  order  to  ensure  complete  coverage  of  the  mesa  step  height.  After  Ni  evaporation,  lift-off  was 
performed  to  remove  the  unwanted  metal  and  resist  using  acetone.  Following  a  thorough  clean 
using  methanol/DI-water,  the  samples  were  again  dehydrated  and  spin-coated  with  a  thick  resist 
layer.  The  samples  were  patterned  using  a  final  metallization  mask  set.  A  layer  of  Cr  and  Au  was 
evaporated  on  the  electrode  regions.  30  nm/200  mn  of  Cr/Au  were  evaporated  and  liftoff  process 
was  used  to  remove  the  resist  and  unwanted  metal.  Figure  3(c)  shows  SEM  pictures  of  a  fully 
fabricated  wall  device. 

3.4  DC  measurements 

At  room  temperature  only  a  small  number  of  carriers  are  thermally  generated  (as  dark 
current)  for  a  Si  bandgap  of  1 . 1 5  eV.  At  low  bias  voltages  (linear  region  of  operation)  the  slope  of 
the  I-V  dark  current  is  proportional  to  the  device  resistance  that  includes  contributions  of  thermally 
generated  carriers  from  both  the  wall  channels  and  the  metal/semiconductor  contact  regions.  At 
higher  biases  the  current  saturates  when  all  thermally  generated  carriers  are  collected.  Any  further 
increase  in  the  current  can  be  attributed  to  leakages  across  the  contact  metal-semiconductor  barrier 
and  to  non-linear  generation  of  carriers  across  the  barrier.  [18]  The  back  interpolation  of  this 
leakage  current  to  the  zero  bias  (0  V)  is  a  measure  of  the  saturated  dark  current  (Ids).  Although  the 
photocurrents  are  a  few  orders  of  magnitude  larger  than  the  thermally  generated  dark  currents,  the 
analysis  of  the  photocurrent  (Ips)  IV  function  is  the  same  as  the  dark  current  (Ids)  IV  plots.  For  dc 
response  analysis,  three  sets  of  measurements  were  perfonned.  These  include  (i)  dark  currents  as 
a  function  of  wall  width  thickness,  (ii)  photocurrents  as  a  function  of  wall  width  thickness,  and 
(iii)  spectral  responses  as  a  function  of  wall  thicknesses.  These  results  are  discussed  and  analyzed 
below. 
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3.5  Dark  currents  versus  wall  width  thickness 


To  study  the  carrier  conduction  properties  versus  dimensionally  scaling  down  the  width 
of  the  wall  structures  into  the  nano-regime,  the  samples  were  characterized  in  batches.  Using 
samples  with  wall  widths  of  200  nm,  95  nm,  75  nm,  40  mn,  and  20  nm,  the  room  temperature 


dark  currents  were  measured  with  a  probe 
station  and  digital  I-V  curve  tracer.  As  the 
physical  cross-sectional  area  of  the  wall 
widths  was  reduced  from  200  nm  to  20 
nm,  we  know  from  Ohm’s  law,  the 
resistance  should  increase  linearly  as  a 
function  of  area.  In  other  words  the 
resistivity  in  units  of  Q/cm2  should 
remain  constant. 
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However,  as  can  be  seen  from  the 
Figure  4,  the  resistivity  is  not  constant  but 
drops  significantly  as  the  width  of  the  wall 


Figure  4.  Plot  showing  resistivity  characteristics 
as  a  function  of  down  scaling  the  wall  widths. 


is  reduced  below  95  nm.  This  suggests  that  there  is  an  increase  in  conductivity  as  the  wall  thickness 
decreases  from  95  nm  to  20  nm.  Since  the  number  of  thermally  generated  carriers  is  directly 
proportional  to  the  volume  of  the  active  region,  any  increase  in  the  conductivity,  as  wall  width 
cross-sectional  region  decreases  from  95  nm  to  20  nm,  cannot  be  attributed  to  the  volume  of  the 
semiconductor  material,  but  must  be  the  result  of  a  substantial  increase  in  the  carrier  velocity. 
Confirmation  of  this  hypothesized  mechanism  was  obtained  with  the  use  of  transient  time  analysis 
as  discussed  in  section  G.  Further  verification  of  this  hypothesized  mechanism  was  performed 
through  transient  time  analysis  as  discussed  in  section  3.8. 

3.6  Photocurrents  versus  wall  width  thickness 

DC  steady  state  photocurrents  were  measured  using  a  365  nm  wavelength,  1.132  W/cm2 
argon-ion  laser  and  a  633  nm  wavelength,  3.96  W/cm2  HeNe  laser.  The  laser  beam  spot  diameter 
was  less  than  8  pm  and  was  focused  within  the  active  region  of  the  electrode  spacing  covering 
several  wall  structures.  By  using  365  nm  and  633  nm  wavelengths,  a  more  complete  insight  into 
absorption  and  carrier  transport  as  a  function  of  wall  thickness  can  be  achieved.  At  365  nm, 
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absorption  occurs  within  the  top  first  10  nm  of  the  Si  wall  structures  with  heights  of  1500  nm.  For 
633  nm  the  total  photon  absorption  extends  through  the  entire  wall  height.  Figure  5(a)  and  5(b) 
show  the  conductivity  versus  wall  thickness  profiles. 


Wall  Width  (nm) 


0  50  100  150  200  250  300  350 

Wall  Width  (nm) 


Figure  5.  Plots  of  photoconductivity  characteristics  as  a  function  of 
down  scaling  the  wall  widths  (a)  for  l=365nm;  (b)  for  l=633nm 


As  can  be  noted  from  the  figures,  a  peak  in  the  conductivity  occurs  around  the  40  nm  wall 
physical  width  of  the  samples  followed  by  a  decrease  in  the  25  nm  wall  width  sample.  The 
significance  of  this  can  be  discussed  through  the  effects  of  strain  on  the  structure  on  the  mobility 
of  the  carriers  as  the  dimensions  are  reduced  in  section  3.8. 


3.7  Transient  time  response  measurements  and  analysis 

The  schematic  of  the  pulsed  carrier  transport  experiment  is  shown  in  Figure  6. 

Pulsed  Laser  Illumination 
+2.5V  hv>E" 


54  I  Si  Nano-walls  si  v 


54  {  Si  Nano-walls 

Gap=  8pm 
Si02  (BOX  Layer) 


Figure  6.  Schematic  configuration  of  a  wall  structured  MSM 
device  used  for  carrier  time  response  measurements 
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When  a  narrow  pulse  of  light  strikes  the  wall  structured  active  region  of  the  device  near 
the  left  electrode  as  shown  in  Figure  6,  equal  number  of  electrons  and  holes  are  generated  and  are 
then  subjected  to  diffusion  and  drift  forces  in  a  presence  of  an  electric  field.  Based  on  the 
experimental  configuration  the  electrons  will  be  rapidly  collected  near  the  positively  biased 
electrode  and  the  holes  will  have  to  travel  the  entire  active  region  to  the  negatively  biased 
electrode.  From  the  measured  time  response  signal  profile  at  the  opposite  electrode,  the  hole 
transient  time  limited  carrier  velocity  can  be  determined,  provided  the  carrier  lifetime  is  greater 
than  the  total  transit  time.  If  the  optical  pulse  of  light  strikes  near  the  opposite  electrode,  the  holes 
will  be  rapidly  collected  and  the  electrons  would  have  to  transit  through  the  active  region,  thus  the 
measured  signal  at  the  opposite  electrode  would  be  electron  transit  time  limited. 

The  pulsed  response  measurements  were  taken  using  150-fs  duration  excitation  at  X  =  400 
nm  from  a  cw  mode-locked  TkAbCb  laser  (doubled  for  the  short  wavelength,  0.2  mW  average 
power  at  a  77  MHz  repetition  rate).  The  wall  structured  MSM  devices  were  probe  tested  using  an 
18  GHz  probe  and  a  high-speed  digital  sampling  oscilloscope  with  an  approximately  1  ps 
resolution  capability.  The  laser  spot  size  was  1  pm  in  diameter  and  the  electrode  gaps  were  8  pm. 
Normal  incidence  was  used  for  the  experiment.  The  time  response  measurements  were  taken  for 
low  electric  field  strengths  3  x  103V/cm,  (2.5  V  across  8  pm  gap)  thus  avoiding  velocity  saturation. 

Before  the  experimental  data  and  analysis  is  provided  it  is  useful  to  review  the  three 
primary  factors  that  can  impact  the  carrier  transport  through  a  semiconductor  region.  These  factors 
are: 

•  Field  dependent  velocity  of  carriers  through  the  active  region.  At  high  E-fields,  the  velocities  of 
both  electrons  and  holes  in  Si  saturate  at  about  1  x  107  cm/s,  [19]  provided  the  field  within  the 
electrodes  exceeds  the  saturation  value  for  most  of  its  length,  we  can  assume  that  the  carriers  move 
with  an  average  velocity  drift.  Velocity  saturation  is  not  an  issue  in  our  experiment  since  the 
applied  field  is  much  lower  than  what  is  required  for  saturation. 

•  Diffusion  of  carriers  in  the  active  region.  The  time  it  takes  for  carriers  to  diffuse  a  distance  d  is 
idiff  =  d2/2D  where  D  is  the  carrier  diffusion  coefficient.  The  diffusion  of  carriers  becomes  a  two 
dimensional  process  as  the  thickness  of  the  Si  wall-structures  is  reduced  and  carriers  are  physically 
constricted  in  movement  by  the  Si/SiO2  interfaces  from  all  sides. 

•  Junction  and  parasitic  capacitance  effects.  A  metal-semiconductor  junction  under  reverse  bias 
exhibits  a  voltage-dependent  capacitance  caused  by  the  variation  in  stored  charge  at  the  junction 
represented  by  the  relation  Cj=A !2(2e SsNd) 1  a V 1 12 ,  where  the  parameters  have  their  usual  meaning. 
This  capacitance  is  usually  quite  small  for  MSM  device  structures  as  a  result  of  their  planar 
electrode  design.  There  are  also  parasitic  circuit  capacitances  associated  with  the  probing  and 
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cabling  that  usually  dominate  the  electrical  response  as  well  as  the  limiting  response  of  the 
electronics.  For  this  study  all  film  devices  have  an  identical  Xdrcuit  limitation. 

Figure  6  represents  the  schematic  with  bias  polarity  of  our  experiment  in  which  the  left 
electrode  polarity  is  positive  and  the  right  electrode  is  ground.  With  this  bias  configuration  once  a 
pulse  of  light  with  a  spot  size  <1  pm,  as  in  the  case  of  our  experiment,  strikes  within  the  active 
region,  the  holes  travel  towards  the  right  electrode  and  the  electrons  travel  in  the  opposite  direction 
towards  the  left  electrode. 

Figure  7  shows  the  experimental  results  of  the  time  response  measurements  for  -200  mn, 
-95  nm,  ~40  mn  and  -20  mn  thick  wall  devices  for  both  electron  and  hole  dominated  signals. 
From  a  first  pass,  as  can  be  seen  from  these  plots,  as  the  thickness  of  the  film  decreased,  the  time 
response  signal  decays  faster.  In  particular  in  the  case  of  the  ~40  nm  and  ~20  nm  thick  walls  the 
signal  decays  over  an  order  of  magnitude  faster  over  the  -200  nm  sample  for  both  electrons  and 
holes.  The  rise  time  is  of  the  signals  is  important  since  it  gives  direct  insight  to  the  carrier  mobility. 
The  rise  time  (td)  is  defined  as  is  the  time-lapse  from  moment  when  the  pulse  of  light  strikes  one 
end  of  the  active  region  of  the  MSM,  near  one  electrode,  and  the  moment  when  the  photo-generated 
carrier  signal  is  detected  at  the  opposite  electrode.  From  the  rise  time  data  we  can  now  determine 
the  carrier  mobilities  as  a  function  of  thickness  as  follows: 

From  the  experimental  time  response  measurements  we  can  calculate  the  average  carrier 
velocities  by  applying  the  given  relation,  V(Carrier- Velocity)  —  (Electrode  gap )/td  (cm/s)  where  td  is  the 
average  time  it  takes  for  the  pulsed  carrier  signal  to  cross  the  electrode  gap  distance.  The  pulse 
travels  in  the  presence  of  a  field  and  expands  from  its  originating  point  due  to  diffusion.  In  this 
case  we  are  ignoring  the  RC  time  delay  that  the  pulsed  signal  experiences  once  it  reaches  the  edge 
of  the  depletion  region  near  the  electrodes  since  the  widths  of  the  depletion  regions  are  very  small 
in  the  submicron  range  compared  to  the  electrode  gap  which  is  ~8pm  in  length. 
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Figure  7.  Measured  time  response  signals  of  200nm  wall  (row-1), 
95nm  wall  (row-2),  40nm  wall  (row-3)  and  20nm  wall  (row-4) 


By  definition  the  average  carrier  mobility  can  be  written  as, 


jU avg  —  V(Avg  Carrier- Velocity)/[VTiias/(EleCtrode)] 


where  VW  is  the  external  bias  applied  to  the  electrodes. 


Figure  8  shows  a  plot  of  average  field  dependent  electron  and  hole  limited  mobility  values 
using  experimental  values  of  rise  time,  td,  and  the  above  expression  as  a  function  of  film  thickness. 
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Figure  8.  Carrier  mobility  values  calculated  from  direct  measure  of 
rise  time  values  as  a  function  of  wall  thickness. 


We  know  that  the  carrier  transport  of  electrons  and  holes  in  the  thickest  wall  sample  (-200 
nm)  is  essentially  similar  to  the  transport  properties  in  bulk  silicon.  However  we  observe  a 
considerable  increase  in  low  filed  dependent  mobility  values  below  -75  nm  wall  thicknesses. 
Recall  the  fact  that  we  actually  have  a  much  narrower  effective  cross-sectional  regions  form  which 
carriers  propagate  due  to  the  repulsive  nature  of  the  boundary  at  the  Si/SiC>2  interface,  and  the 
carrier  profile  tends  to  peak  a  certain  distance  away  from  the  interface  close  to  the  center  of  the 
wall  structures.  Therefore  we  believe  that  quantum  confinement  effects  play  a  key  role  in  the 
carrier  mobilities  of  both  electrons  and  holes.  A  detailed  theoretical  model  is  described  in  the  next 
section  which  tends  to  explain  our  experimental  results. 


3.8  Theoretical  Modeling  using  strain  effects  to  explain  rise  in  electron  and  hole 
mobility 

If  we  consider  a  total  valence-band  hole  concentration  nv  then,  the  light-hole  (nLH)  and 
the  heavy-hole  (nHH)  concentration  will  satisfy  the  charge-conservation  relation  nLH  +  nHH  —  nv 
,  where 


9r9s 


[l  +  exv  *  2  gr  (^ff2  exp  , 
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where  the  subscript  er  takes  HH  or  LH  and  the  upper  (lower)  sign  corresponds  to  HH  (LH)  state. 
In  the  above  expressions,  the  approximations  are  made  for  high  temperatures,  V  is  the  volume  of 
the  silicon  film,  T  is  the  system  temperature,  the  zero  energy  is  chosen  at  the  middle  point  between 
the  split  pair  of  light-hole  and  heavy-hole  bands,  k  is  the  three  dimensional  wave  vector  of  carriers, 
gr  —  2  (not  6  due  to  strain  effect)  is  the  F-vallcy  degeneracy  for  holes  and  gs  —  2  is  the  spin 
degeneracy  for  both  light-holes  and  heavy-holes.  In  addition,  uv,  which  depends  on  both  T  and 
nv,  is  the  chemical  potential  to  be  determined  for  valence  bands,  E%h  —  h2k2 /2m*HH  is  the  kinetic 

energy  of  heavy  holes  and  EfcH  —  — —  is  the  kinetic  energy  of  light  holes,  where  m*HH  — 

2mLH 

0.49m0  and  m*LH  —  0.16 m0  (m0 is  the  free-electron  mass)  are  the  effective  masses  for  heavy  holes 
and  light  holes,  respectively.  Additionally,  A Egtr  introduced  in  the  above  expressions  stands  for 
the  half  of  the  valence-band  splitting  due  to  the  existence  of  strain. 

From  the  above  two  equations  and  nLH  +  nHH  =  nv,  we  obtain  nLH/nv  =  [1  + 
y3/2  expilgAEsfy/kgT)]-1  and  nHH/nv  —  1  —  nLH/nv,  where  y  =  ( [m*HH/m*LH )  >  1  .  For 
biaxial  and  shear  strains,  [1,2]  we  have  the  valence-band  splitting,  given  by  A Egtr  = 

2  2  1  /2 

i  {(b  /2)  £yy)  +  (^yy  ^zz)  +  C^zz  ^xx)  J  T  d  \^xy  "F  £ yz  "F  ^xz\ |  ■>  where  the 

upper  sign  is  for  the  compressive  strain  while  the  lower  sign  for  the  tensile  strain  in  the  direction 

perpendicular  to  the  interface  of  silicon  and  silicon-dioxide  materials,  b  and  d  are  the  optical 

deformation  potentials,  and  *  represents  the  strain  tensor  in  the  three  dimensional  space  with  j, 

j’  =  x,  y,  and  z,  the  diagonal  matrix  elements  ejj  are  associated  with  biaxial  strain,  and  the  off- 

diagonal  matrix  elements  with  j  ^  / ’  correspond  to  contributions  from  the  shear  strain.  For 

silicon  crystals,  we  have  b=  -2.33  eV  and  d=  -4.75  eV. 

The  values  of  g  can  be  scaled  as  g  =  (l/L){min(L,  2 Ds)} ,  where  L  is  the  film  thickness 
and  Ds  is  the  strain  buffing  depth  due  to  lattice  mismatch  between  embedded  Si  crystal  and 
surrounding  amorphous  SiCk  material  at  their  interface,  and  L  —  2 Ds  >  0  represents  the  film 
effective  thickness  free  of  localized  trapping  centers  [19]. 

If  we  choose  the  z  direction  as  the  direction  perpendicular  to  the  interface  for  biaxial  strain 
we  simply  get  exx  =  eyy  =  e„ ,  ezz  =  e±,  and  =  0  for  i  A  j,  [19]  wheree,,  =  ( alsi/asi  -  l), 
=  ( a±,si/asi  ~  l).  Moreover,  the  perpendicular  lattice  constant  a±si  is  related  to  the  parallel 
lattice  constant  alsi  =  asi02  by  a±si  =  asi[l  -  (2c12/c11)(asi02/asi  -  l)],  where  clt  = 
16.75  x  1010  N/m2,  and  c12  =  6.5  x  1010N/m2  are  the  elastic  constants  of  silicon.  For  silicon 
and  silicon-dioxide,  we  have  asic>2  =  (2  x  4.914  +  5.405)/3  =  5.078  A  and  asi  =  5.431  A  for 
amorphous  silicon-dioxide  materials.  Therefore,  we  obtain  a±si/aSj  =  1.050.  This  leads  to  f||  = 
—0.065  (compressive),  =  0.05  (tensile),  and  2 ey  +  e±  —  —0.08. 
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The  total  mobility  )tv  for  holes  can  be  expressed  as 


'nw') 

|  eTLH 

.  nv  ) 

m0 

m0  ^  /  m0  \1  / nHH\  erHH  wip  ^  /  m0  \1  _  / nLH\  ezLH  +  / nHH\ 

mLH  'jrlLH' J  V  nv  /  m0  L^hh  VmhH'J  V  nv  /  \  nv  J 


ezHH 


m 


HH 


lLH 


‘HH' 


(l+y3/2) 


fezLii 

6THh\ 

exp(2?7A£-ftr//cBr)-l 

^mLH 

™hh' 

(l+y3/2)2 

(0) 

kv 


1  - 


??(2A£'ftr/fcBr)y2(y-i) 


(l+y1/2)(l+y3/2) 


+  0 


©  ■ 


fe.,  -  l)  oc  r;  oc  1/L,  where  /4,0)  =  (eTw/m*LH)(l  +  T1/2)(l  +  73/2)  is  the  valence 

band  mobility  for  r\  0,  l/fv  =  (1/2)(1/tlh  +  l/rHH)  (withrLH  *  rHH  «  fv), \AE”tr\  «  kBT 
is  assumed,  the  changes  in  the  hole  effective  masses  by  strain  have  been  neglected,  rLH  and  rHH 
are  the  scattering  times  for  light  and  heavy  holes,  respectively.  It  is  clear  that  fiv  increases  with 
1/L  for  the  tensile  strain  (A Egtr  <  0)  in  the  direction  perpendicular  to  the  interface  of  silicon  and 
silicon-dioxide  materials,  as  observed  by  us  in  Figure  9. 


Figure  9.  Theoretical  modeling  for  electron  (left  panel)  and  hole  (right  panel)  mobilities  ias  functions  of 
wall  thinckness  L  with  a=  1.0  (red  solid  curves)  and  1.5  (black  dashed  curves)  and  their  comparisons 

with  experimental  data  (black  dots  in  both  panels 

d  electron  concentration  nc,  the  electron  chemical  potential  uc,  which  depends  on  both  T  and  nc, 
is  decided  from 


Approved  for  public  release;  distribution  is  unlimited. 

16 


where  the  high-temperature  approximation  is  made  in  the  above  expression,  Eq  —  £q(T )  +  77' A Eq 
is  the  bandgap  energy  of  strained  silicon  crystals,  which  depends  on  T  and  the  hydrostatic  part  of 
the  strain,  £q  stands  for  the  bandgap  energy  of  unstrained  silicon  crystals,  gXL  —  2  (not  6  due  to 
strain  effect)  represents  the  X  (in  <100>  direction)  and  L  (in  <1 1 1>  direction)  valley  degeneracy 
for  electroncs  at  the  two  minima  of  conduction  band,  E^—h  k2 /2m jE  is  the  kinetic  energy  of 
electrons  and  m ^  is  the  transverse  effective  mass  of  conduction-band  electrons  with  m*x  —  0.19m0 
and  m*L  =  0.1m0.  The  T  dependence  of  £q(T)  (based  on  the  Bose-Einstein  phonon  model)  is 
given  by  £q(T)  =  £q(0)  —  2aB@B[coth(QB/2T )  —  1],  where  aB  =  2.82  x  10-4eV/K  is  a 
coupling  constant,  kB&B  is  a  typical  phonon  energy  with  @B  —  351  K,  £q(T)  —  1.12  eV  and 
£q(T)  —  2.4  eV  at  T  —  300  K  for  the  X  and  L  valleys.  Moreover,  the  strain  part  [2]  of  the  bandgap 
energy  A Eq  is  calculated  as  A Eq  —  E/pTrie)  +  E/pe^  •  €*•  e^  +  aT r (e)  ,  where  EpL^  and  e/X,L'> 
are  the  deformation  potentials  of  the  conduction  band  for  an  indirect-gap  silicon  crystal  (Ed  = 
1.1  eV,  E(p  —  10.5  eV  for  the  X  valley  and  E(p  =  —  7.0  eV,  Ep*  —  18.0  eV  for  the  L  valley), 
a  —  2.1  eV  is  the  difference  of  the  deformation  potentials  of  conduction  and  valence  bands  at  two 
different  valleys  due  to  hydrostatic  component  of  the  strain  for  the  silicon  crystal,  and  e ^  is  the  unit 

vector  pointing  to  the  specific  X  or  L  valley.  It  is  clear  from  the  above  equation  that  A Eq  <  0  for 
the  tensile  strain  and  ^  =  X  or  L  . 


The  change  in  the  bandgap  energy  by  strain  also  affects  the  effective  mass  of  conduction 
band,  given  by  [20]  also  see  [21-33]. 


A 


£p(2e  ii+g-OV 

4(7’)+ A0/3 


2  + 


4(r)+A0/3 


+  0[(2f||  +  £±)2]  , 


where  we  have  neglected  the  shear  strain  and  assumed  a  weak  strain  with  \2e^  +  £±|  «  1  ,  A0  = 
44  meV  is  the  spin-orbit  splitting  and  EP  —  21.6  eV  is  the  Kane  energy  parameter. 


The  total  mobility  gc  of  conduction-band  electrons  is  obtained  as 


+ 


£f![2!£  +  a 

m0  Ym*x 


+ 


_ fix Yi  4.  ( exA  Ep(2e ll+e-i) 

m*x  ^  \m0)  £q(T)+A0/3 

=  ,<<0)  [l  -  fj'Cl  +  «)  g) 


m0J  £*(T)+ Ao/3  [ 

Ep(2  ell+e±) 


3a+3=-m+=.w 

2  _|_  ia+5~d  +“u 


£§(T)+A0/3 


4(T)+A0/3 


n  3a+34X)+£(X)l 
2  + - 4 - 


4(T)+A0/3 


+  <3[(2£||  +  £±)2] 

|  +  0[(2f||  +  e±)2] , 
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i.e.,  (\ic/\Pc  —  1)  oc  77'  +  0(t]'2)  oc  1/L  +  0(1/L2)  (using  rj' ~  r]),  where  we  assume  = 
r°(m0m|)a  with  a  labelling  the  mass  dependence  in  the  scattering  rate,  fi(c0>  «  er x/m-x  is  the 
conduction-band  mobility  for  rj'  ->  0,txl  represents  the  scattering  times  of  conduction-band 
electrons  at  two  different  valleys  and  the  high-energy  L  valley  has  been  assumed  depopulated.  It 
is  clear  that  the  electron  mobility  increases  with  1/L2  in  our  system  with  (2e\\  +  e±)  =  —0.08,  as 
observed  by  us  in  Figure  9. 
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4  RESULTS  AND  DISCUSSION 


Table  1  and  2  below  show  the  results  of  the  theoretical  calculations.  As  was  previously 
shown  in  Figure  9,  the  comparisons  between  experiment  and  theory  are  within  10%  which  is  a 
very  good. 

Table  1.  Model  parameters  used  in  calculating  mobility  of  electrons  in  strained  Si 
WALL  Structures. 


i/max  2 

c  (cm  /V  ■  s ) 

fJc0)  (cmr/V  ■  s) 

A  (nm) 

2  Dc  (nm) 

5500 

806 

15 

42 

Table  2.  Model  parameters  used  in  calculating  mobility  of  holes  in  strained  Si 
wall  structures. 


U  max  2 

v  (cm  /V  ■  s ) 

fjy  ])  (i cm2/V  ■  s) 

/\v  (nm) 

2 Dv  (nm) 

3000 

100 

143 

42 

As  can  be  noticed  from  both  the  experiment  and  theory,  the  carrier  concentration  includes 
both  the  doping  and  photo-excitation  contributions.  If  the  sample  is  undoped,  we  can  simply 
neglect  the  impurity  scattering  and  have  nc  —  nv.  The  optical-phonon  scattering  and  the  inter¬ 
valley  scattering  are  only  important  at  high  temperatures,  while  the  acoustic-phonon  scattering 
becomes  more  significant  at  low  temperatures.  The  surface-roughness  scattering,  on  the  other 
hand,  is  largely  independent  of  temperature.  In  the  narrowest  wall  devices,  a  considerable  increase 
in  conductivity  was  observed  as  a  result  of  higher  carrier  mobilities  due  to  lateral  constriction.  The 
strain  effects,  which  include  the  reversal  splitting  of  light-  and  heavy-  hole  bands  as  well  as  the 
decrease  of  conduction-band  effective  mass  by  reduced  Si  bandgap  energy,  are  formulated  in  our 
microscopic  model  for  explaining  the  experimentally  observed  enhancements  in  both  conduction- 
and  valence-band  mobilities  with  reduced  Si  wall  thickness.  Specifically,  the  enhancements  of  the 
valence-band  and  conduction-band  mobilities  are  found  to  be  associated  with  different  aspects  of 
physical  mechanisms.  The  role  of  the  biaxial  strain  buffering  depth  is  elucidated  and  its  importance 
to  the  scaling  relations  of  wall-thickness  is  reproduced  theoretically,  i.e.,  1/L2  for  electrons  and  1/L 
for  holes. 
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5  CONCLUSION 


The  semiconductor  processing,  fabrication  and  the  resulting  carrier  transport  character¬ 
istics  of  MSM  devices  fabricated  as  wall  like  structures  in  silicon  on  insulator  technology  were 
reported.  MSM  device  dark  current,  DC  photocurrents,  and  the  time  response  of  car-rier  transport 
were  investigated.  The  resulting  conducting  channels  were  actually  smaller  than  their  physical 
dimensions,  a  result  of  depletion  of  carrier  near  the  interfaces.  As  the  physical  channel  widths  were 
reduced  by  oxidation,  strain  was  produced  near  the  interface  and  strained  lattice  became  a 
significant  portion  of  the  conducting  channel.  The  increase  in  mobilities  for  both  holes  and 
electrons  stemming  from  the  strained  silicon  resulted  in  a  dramatic  increase  in  carrier  mobility  for 
both  electrons  and  holes  as  the  physical  channel  width  was  reduced  from  200  nm  to  20  nm.  The 
theoretical  model  incorporating  the  effects  of  strain  present  in  these  nanoscale  MSM  devices 
compared  favorably  with  experimental  results,  showing  that  hole  mobilities  increased  with 
decreasing  L.  Additionally,  if  these  elec-tron  and  hole  mobilities  can  be  retained  with  the 
application  of  gate  electrodes,  then  this  technique  may  yield  a  much  simpler  path  towards  high 
performance  CMOS,  both  n-channel  and  p-channel,  than  current  techniques  for  either  planer  ultra- 
thin  body  FETs  or  FinFETs. 
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LIST  OF  ACRONYMS 


BOE 

Buffered  oxide  etch 

CMOS 

Complementary  Metal  Oxide  Semiconductor 

DC 

Direct  current 

FET 

Field  Effect  Transistor 

FinFET 

Fin  shaped  Field  Effect  Transistor 

FM 

Fresnal  Mirror 

HH 

Heavy  hole 

Ids 

Dark  current 

IL 

Interferometric  Lithorgraphy 

Ips 

Photocurrent 

ITRS 

International  Technology  Roadmap  for  Semiconductors 

LH 

Light  hole 

MSM 

Mobile  Service  Manager 

MSM 

Metal  Si/nanowall  Metal 

nFET 

n-type  FET 

pFET 

p-type  FET 

Qf 

Fixed  Charge 

Qa 

Interface  Trapped  Charge 

QoT 

Oxide  Trapped  Charge 

RIE 

Reactive  ion-etching 

SEM 

Scanning  Electron  Microscope 
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SGOI 

SOI 

sSOI 

VDD 

VLSI 


SiGe  on-insulator 
Silicon-on-Insulator 
Strained  Silicon-on-Insulator 
Voltage  Drain  Supply 
Very  Large  Seal  Integration 
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